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• ABSTRACT 

II ■ Context. A stellar wind passing through the reverse shock is deflected into the astrospheric tail and leaves the stellar system either as a sub- 

^ I Alfvenic or as a super-Alfvenic tail flow. An example is our own heliosphere and its heliotail. 

}—i . Aims. We present an analytical method of calculating stationary, incompressible, and field-aligned plasma flows in the astrotail of a star. We 

C/J present a recipe for constructing an astrosphere with the help of only a few governing parameters, like the inner Alfven Mach number and the 

^ ] outer Alfven Mach number, the magnetic field strength within and outside the stellar wind cavity, and the distribution of singular points (neutral 

points) of the magnetic field within these flows. 

^_^ Methods. Within the framework of a one-fluid approximation, it is possible to obtain solutions of the governing MHD equations for stationary 

^ flows from corresponding static MHD equilibria, by using noncanonical mappings of the canonical variables. The canonical variables are the 

C^ Euler potentials of the magnetic field of magnetohydrostatic equilibria. Thus we start from static equilibria determined by the distribution of 

^^ magnetic neutral points, and assume that the Alfven Mach number for the corresponding stationary equilibria is finite. 

'/^ , Results. The topological structure, i.e. the distribution of magnetic neutral points, determines the geometrical structure of the interstellar gas - 

^'i • stellar wind interface. Additional boundary conditions like the outer magnetic field and the jump of the magnetic field across the astropause al- 

fT^ ' low determination of the noncanonical transformations. This delivers the strength of the magnetic field at every point in the astrotail/astrosheath 

f^ region beyond the reverse shock. 

OQ Conclusions. The mathematical technique for describing such a scenario is applied to astrospheres in general, but is also relevant for the helio- 

^~~^ , sphere. It shows the restrictions of the outer and the inner magnetic field strength in comparison with the corresponding Alfven Mach numbers 

S>. ' in the case of subalfvenic flows. 
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1. Introduction function, and the method is restricted to 2D fields and flows. 

We apply a method that is not restricted to one flux function, 



1.1. The scenario 



but works with two flux functions of the magnetic fields. 



A wide range of literature is concerned with calculating sta- 
tionary MHD flows for stellar magnetospheres, jets, stellar Here we show that, under certain reasonable assumptions, 

winds, and laboratory or general plasma configurations, see it is possible to use a powerful transformation method for sys- 

e.g. Chandrasekhar (fT956] l. Tsinganos (TWPt, Lovelace et al. tematic modelling of the stellar wind region far away from the 

(fT986] l. and Goedbloed & Lifschitz ( 1997). These authors use star itself. We apply this method to the special scenario of a 

one flux function to represent two components of the magnetic stellar wind-interstellar medium (ISM) counterflow configura- 

field and get a Grad-Shafranov type equation. This is a single tion (see Fig.fflfor the special case of the heliosphere). From 

nonlinear, partial diff'erential equation for this magnetic flux this scenario, it is possible to estimate the pressure of the mag- 

netic field, plasma pressure, and ram pressure, which are dy- 

Send offprint requests to: D.H. Nickeler namically important for the ISM and therefore of high interest 

Correspondence to: nickeler@asu.cas.cz to astronomy (see e.g. Frisch[T993}. 
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Fig. 1. Sketch of the heliosphere shown as a special example of 
an astrosphere. 

Beyond the region of the reverse shock ('Heliospheric 
Shock' in Fig.[TJ, the plasma of the stellar wind is decelerated. 
The magnetic field increases at the shock, so that a sub- or a 
super- Alfvenic plasma flow exists in the downstream direction. 
A contact surface forms between the two different flows. This 
separatrix is called the astropause. Scherer et al. (,1994 1 and 
Fahr et al. (I1993I I showed that the bulk flow in the downwind 
direction can be assumed to be incompressible for small Mach 
numbers. This holds even more for field aligned flows, since 
the field lines act as quasi-isothermals. As the decelerated stel- 
lar wind has to adapt to the conditions of an outer magnetized 
Very Local InterStellar Medium (VLISM), a tangential discon- 
tinuity forms. This is called the astropause (AP, indicated by 
'Heliopause' in Fig.[T]), which stretches out in the downwind 
direction, so that the whole structure has a tail-like shape (see 
Fig.O. Between the termination shock (TS, the heliospheric 
shock in Fig.lTJ and the AP an inner astrosheath region extends 
into the astrotail, similar to the Earth's magnetotail. In the mag- 
netotail, the MHD quantities depend mainly on the direction 
perpendicular to the tail axis (e.g. Schindler ll972) . In view of 
the evident similarities, we here apply a similar description to 
astrotails. We show that there is a strong correlation between 
the flow, especially its Alfven Mach number inside and outside 
the astropause region, and the current density in the vicinity of 
the astropause. This correlation also sets restrictions on the re- 
lation between outer and inner values of the magnetic field and 
the Alfven Mach numbers. 

1.2. Model considerations and existing models 

Astrospheric models can be divided into two kinds, analyti- 
cal and numerical, but there is a gap between these two. In 
most cases, the analytical models are only hydrodynamical 
or purely magnetic models (e.g. Parker[T96T]. Parker's article 
only treats purely hydrodynamic models of the subsonic coun- 
terflow or of an unmagnetized stellar wind that blows into a 
magnetohydrostatic interstellar environment (Parker[T96T]l. His 



models connect parameters, like the pressure of the magnetic 
field and the thermal pressure , with the shape of the model- 
astrospheres. Another analytical model was the first super- 
Alfvenic MHD model suggested for the SW-ISM interaction, 
which is based on the thin-layer (hypersonic) approximation 
(Baranov & KrasnobaevlWT"). In contrast to the model of 
Parker or Baranov and Krasnobaev the geometry of our model, 
presented here in this article, does not explicitly depend on the 
Alfven Mach number or the usual Mach number of the flow. 
We propose a different point of view in this article by emphasiz- 
ing the connection between topological aspects of the magnetic 
field structure and the geometrical shape of the astropause. 

Other analytical models have been calculated for pressure 
equilibrium (Newtonian approximation by Fahr et al. 1988) or 
for those configurations where the plasma cavity is a finite el- 
lipsoid, and the plasma has to leave the astrosphere by diffu- 
sion (Neutsch & Fahi 11983l l. Other authors use the set of MHD 
equations, but prefer solving the ideal Ohm's law and neglect- 
ing the Lorentz force in the Euler equation (see the kinematical 
approach by Suess & Nemev lT992llT990l Nerney et al. [T99Tl . 
Up to now no analytical and exact solutions of the MHD equa- 
tions exist for this scenario. This motivated us to consider mod- 
els that do not depend on mathematical approximations, but 
where (additional) physical approximations are taken into ac- 
count. For example, Imai ( 119601 1 analysed field-aligned flows 
and calculated approximative solutions. However, we take into 
account that a flow that has passed a shock is likely to develop 
a sub-sonic/sub-Alfvenic flow with a negligible compressibil- 
ity along the field- and streamlines. The argument in Scherer et 
al. ( II 9941 1 and Fahr et al. (11993') is that the low Mach number 
does not provide high compressibility rates. Another argument 
for 'incompressibility' is that in a tail that is symmetric with 
respect to the tail axis, the far-away field- and streamlines are 
nearly one-dimensional. This implies a one-dimensional de- 
pendence of the physical values perpendicular to this axis, see 
Schindler ( 1972) , so that the density is approximately constant 
on field lines. 

Webb et al. ( 119941) analysed two-dimensional MHD flows 
and discussed the properties of transsonic flows. We focus on 
the relation between the Alfven Mach number Ma and the elec- 
tric current density, which can be given by the jump of the two 
magnetic fields across the boundary of the astropause of, in 
principle, 3D configurations. We also focus on topological and 
geometrical questions with respect to the boundary between 
two magnetized flows. 

Simulation results may also be reliable, although they can 
deliver unphysical results, even when stable algorithms are 
used. For example, Linde et al. (119981 ) discovered magnetic dif- 
fusion in their simulation domain, although they used an ideal 
MHD code. Evidently, numerical magnetic reconnection is tak- 
ing place in some region of their domain of calculation. 

Numerical reconnection cannot take place in our treatment 
since we calculate exact and analytical solutions of the ideal 
MHD equations with finite width of the astropause current 
sheet. We describe smooth flows without shocks, i.e. without 
non-tangential discontinuities, only. In our approach, we have 
restricted the analysis to tangential discontinuities. The possi- 
bility of Alfvenic discontinuities is discussed by Smith (120011 ) 



D.H. Nickeler et al.: MHD flows in astrotails 



for the heliospheric current sheet. Thus, we close the gap be- 
tween the older analytical, but simplified, treatments of the as- 
trophysical counterflow problem and the sophisticated numer- 
ical models. We find principles for these counterflow scenar- 
ios that should also hold for simulations having fewer physical 
simplifications. Our aim is to calculate analytical and exact so- 
lutions of the stationary nonlinear MHD equations for stellar 
wind tails, extending the approximative analytical models dis- 
cussed above. 

The paper is structured as follows. In Sect. 2, we present the 
equations to be solved for our scenario and the special method 
used to solve them. In Sect. 3, we show how this method works 
in the two-dimensional case, which is a simplification of our 
model, focussing on the method and the basic principles of 
MHD counterflows. In Sect. 4, we show how to construct the 
flow and its stream lines, i.e. the magnetic field structure. This 
gives us the pattern of such MHD flows. We use the simplest 
magnetohydrostatic equilibria, viz. potential fields. The reason 
for it will be given in Sect. 4. Section 5 discusses the depen- 
dence on the symmetry and boundary conditions for the asymp- 
totical one-dimensional case, which is interesting for tail-like 
structures of the magnetic field. In addition, we present some 
two-dimensional tail models. Discussion and conclusions are 
given in Sect. 6. 

2. Stationary states in incompressible and ideal 
MHD 

The set of equations that must be solved to get incompress- 
ible ideal MHD flows consists of the mass continuity equation 
([T]l, the Euler or momentum equation with isotropic pressure P 
(|2|l, the induction equation including the ideal Ohm's law (|3]l. 
Ampere's law (|4]l, the initial condition for the magnetic field 
(|5]l, and the condition for incompressibility (l6}: 



V-(pv) = 0, 

p (v ■ V) V = jxB 

Vx(vxB) = 0, 

VxB = jjoj, 

V B = 0, 

Vv = 0. 



VP, 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 



Due to the incompressibility, the mass continuity equation can 
be written as v ■ Vp - 0, so that the density is constant on 
streamlines. If we now introduce the auxilliary flow vector w := 
■\/p\ and the Bernoulli pressure H :- P + \w^, we can write 
the above equations as 



Vw = 0, 



(7) 



Vx| — wxB = 0, 
VB = 0. 



Vn = — (V X B) X B - (V X w) X vv , (8) 



(9) 
(10) 



Hence, the momentum equation Eq. ([8]) is written such 
that the analogy with magnetohydrostatic equilibria, VP = 
//g' (V X B) X B, is evident. 



2.1. Field-aligned flows 

The stationary equilibria should be constructed such that they 
tend to be stable in order to use them as stationary background 
fields in very turbulent and time-dependent stellar winds. In 
analytical works (e.g. Suess & Nemey[T992] and references 
therein), this problem is often treated kinematically, which 
means that the Lorentz force is ignored. These authors find 
strong amplification of convected magnetic fields in the so- 
called upwind direction, which is the direction from which the 
interstellar medium is flowing towards the star. In this direction 
they identified a cone of 30 degrees (where the star is sitting at 
the top of the cone), in which their kinematical approach is in- 
valid. Such velocity fields with a strong perpendicular compo- 
nent to the magnetic field have a saddle-point structure in lin- 
ear stability analyses (Hameiri ll998t and are, therefore, likely 
to develop ideal MHD instabilities. 

Hameiri found that a variational principle does not lead to a 
stability criterion if velocity and magnetic field are not aligned, 
because the used functional has only stationary points, but has 
no minimum. Hameiri (1998) suggests that the lack of a mini- 
mum is due to the presence of ballooning modes. In fact assum- 
ing the incompressible limit, the equilibrium velocity field has 
to be sub-Alfvenic to ensure the existence of a minimum. Thus, 
to calculate magnetohydrodynamic configurations that should 
'survive' long enough to represent a quasi-stationary state of 
the stellar wind flow, it is necessary to assume field aligned 
flows. This ensures that these configurations can really exist 
in nature, i.e. that they are sufficiently long-lived to be repre- 
sented as stationary MHD flows. This would not be the case 
with strong perpendicular components of the flow with respect 
to the magnetic field, since those would lead to quick ideal 
MHD instabilities. Thus, models where magnetic and flow field 
are not (approximately) aligned cannot exist in nature without 
showing strong time dependency. Since strong time-dependent 
MHD flows, instabilities, and shocks do occur in the corotat- 
ing interaction regions of the solar system, the validity of our 
models is restricted to those regions far outside the termination 
shock. There, incompressible and field-aligned flows are good 
approximations of the real outflows of stars, at least when they 
are sufficiently far away from the stellar surfacel^ 

From this point of view, to lower the risk of instability, it is 
expedient to make the simplifying assumption of field aligned 
flow, i.e. 



±^B, 



-v^ 



(11) 



where Ma is the Alfven Mach number This equation fulfills the 
induction equation (|9]l automatically. The sign on the righthand 
side of Eq. ( fTTT i is to be understood in the framework of the 
transformation method introduced in the next subsection. With 
this assumption, we can skip Eq. ^, and from Eqs. (fTTT i. ^, 



This is also valid, if we focus our view either on the classical (sub- 
sonic unmagnetized flow, radially extrapolated to the origin) Parker 
flow or on the Parker spiral magnetic field, which we use as a geo- 
metrical pattern for calculating flows and magnetic fields in the next 
sections. 
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and ( [Tol l it follows directly that 

w ■ VMa = , and B VMa = . 



(12) 



Therefore, the mass density p and the Alfven Mach number 
are constant on field lines, but they can vary perpendicular to 
them. In conclusion, we have to solve the following system of 
equations: 

B-VMa = 0, (13) 

Vn = — (l - mI)(V X B) X B - ^—^V (l - Ml) ,(14) 

V-B = 0. (15) 

This system determines the unknowns B, Ma and n. In the 
next subsection we present a method of solving the general 
three-dimensional problem given by the system of Eqs. (fTJt- 
( fTsT l by means of a noncanonical transformation. Later on, we 
explicitly calculate two-dimensional equihbria in order to high- 
light the main properties that are important for understanding 
relaxed, magnetized, stellar tail flows. The reduction of this 
system to one magnetic flux function (e.g. Tsinganos 119811 
Goedbloed & Lifschitz lT997] l will be done in Sect.[3] 

2.2. Euler potential representation and noncanonical 
transformations 



In most cases in the literature (e.g. Chandrasekhar [T956] 
Tsinganos' 1981; Lovelace et al. 119861 Goedbloed & 
Lifshitz 1997) , the problem of solving the stationary MHD 
equations is reduced to equations similar to the Grad-Shafranov 
equation (GSE, see e.g. Grad & Rubin' 195 8) by introducing 
two-dimensional flux functions for the magnetic field. Here, 
we give a short introduction to a different method that allows 
us to calculate also fields that could be three dimensional. 

In 1984, Zwingmann showed the similarity between mag- 
netohydrostatic (MHS) equilibria and stationary MHD equi- 
libria with incompressible, field-aligned flows. Later, this the- 
ory was improved by Gebhardt & Kiessling (1992) , and sub- 
sequently used by Petrie & Neukirch (1 19991 1 for modelling 
sunspot magnetic fields with plasma flow. We briefly recapitu- 
late the transformation method used in the cited papers in order 
to facilitate the analysis of our astrospheric model. 

In general, non-ergodic magnetic fields can be represented 
by using Euler potentials (see e.g. Kruskal & Kulsrud 1958 or 
D'haeseleer l 19901 and references therein). The magnetic fields 
of MHS equilibria can also be represented by using Euler po- 
tentials, writing 



Bs^VfxVg, 



(16) 



where the Euler or Clebsch potentials / and g are scalar func- 
tions of X, y, z in general. Here, and in the following, the sub- 
script S will be used to indicate magnetohydrostatic equilib- 
rium quantities. The MHS equations can now be written as 
canonical Hamiltonian equations: 



dPs . ^ 
-^- = Js -V/, 



(17) 



(18) 



with the canonical variables / and g and the arc length s along 
the current js (see Schindler ll9791 l. The MHS field Bs can now 
be mapped to a new field B by performing the transformation. 



/ = f(a,/3) 



g = g(a,/3) 



a = a(f,g) 



/3^/3(f,g) 



(19) 



where the derivatives as well as the inverse mappings and its 
derivatives, are assumed to exist. Then, there is a relationship 
between the old (static) field Bs and the new field B, which can 
be interpreted as a stationary field: 



Bs^VfxVg^ [/, g]a,pVa X V/3 = [/, g],,^ B , 



(20) 



where the Poisson bracket [/, g]afi is the Jacobian of the trans- 
formation Eq. (fT9] l. If the Poisson bracket [f,g]a.0 = 1, it can 
be seen from Eq. ( l20l l that only a change of the potentials but no 
real active transformation has taken place, so that the magnetic 
field has not been changed. Therefore, a and /? would also be 
canonical variables for the field Bs , and the mapping from Bs 
to B would be a canonical transformation, which does not pro- 
duce new physics. However, if the Poisson bracket [/, g]afi + 1, 
then the magnetic field B has to be interpreted in a diff'erent 
way, as it is not possible to identify it as a magnetic field of a 
static equilibrium. 

The similarity between MHS, Eq. (I2T]| below, and MHD, 
Eq. (fT4l i, can be seen by inspecting the original momentum bal- 
ance equation of the MHS field, when we insert Eq. ( fT9] l into 
Vf s = is X Bs , which leads to 



VP., 



1 (V X (V/ X Vg)) X (V/ X Vg) 

[/ §]^ 

"'^ V X (Vff X VyS) X (Va X Vy3) 

IVaxVySP 2 



2;Uo 



(21) 



Consequently, the relation between the Poisson bracket and the 
Alfven Mach number is given by 



0<([/,gU) :=1-M2 

for purely sub-Alfvenic flows, and by 



(22) 



<([/,§]„/:= Ml -1 (23) 

for purely super- Alfvenic flows. 

Setting B = in Eq. (O, which describes a purely in- 
compressible stationary hydrodynamical flow, Gebhardt & 
Kiessling ( 119921 ) noted the similarity between Eq. ^ and 
Eq. (l2n i. This can be seen if the auxiliary flow field w is 
also represented by Euler potentials. Then a mapping from a 
known solution of stationary incompressible hydrodynamics to 
a stationary incompressible, super-Alfvenic field-aligned flow 
is possible, if {\f,g\afi)^ > 1 . It is possible to map a known 
solution of the MHS equations by means of a transformation 
with ([/, §]q.,^)^ < 1 to an incompressible MHD equilibrium 
with a field-aligned sub-Alfvenic flow. Thus for every incom- 
pressible field-aligned sub-Alfvenic flow, it is possible to find a 
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mapping onto a MHS equilibrium. Alternatively, it is also pos- 
sible to take a different look at the problem by focussing on the 
'current-generating' transformation of a given MHS equilib- 
rium, as we are interested in the astropause current sheet. This 
works as follows: from Eqs. (ISTT i and (l23l l it is obvious that for 
a certain transformation, e.g. given by a sub-Alfvenic flow with 
the Jacobian squared ([f,g]a.i3)^ < 1, 



w,„t = sign[[/,g]„^] 



^|^-{[f^8Uf 



-v^ 



B, 



(24) 



there exists a corresponding super-Alfvenic solution for the 
flow fields; 



= sigi 



n[l/'^]<,,/;] 






+ 1 



B. 



This has to obey the following restriction 



< [f,g]it 



M- 



A.super 



1 = 1- ^Isub < 1 



M- 



A,super 



= 2-M- 



A, sub ' 



to be satisfied at every point in space, while the magnetic field 
is the same as in the sub-Alfvenic case: B super - Bsub = B. 

Therefore, it is not guaranteed for all Poisson brackets, 
i.e. transformations, that sub-Alfvenic solutions exist, but for 
a given sub-Alfvenic solution ([f,g]afi)^ < 1 a corresponding 



super-Alfvenic MHD flow exists with 1 < Afj ^^ < 2. For 
these flows, the magnetic field of the underlying MHS equilib- 
rium magnetic field will be amplified: 



\B\ = Bs 



Bs 



'\j-"^A,super 



>\Bs 



(28) 



It is also necessary for the super-Alfvenic case that the ther- 
mal or plasma pressure Tlsuper is 'inverted' to regain the simi- 
larity between Eqs. ( fT4l i and (ISTT i. 



changed to 

VPs =V [-Tlsuper) 



(29) 



(30) 



where Hsub is the sub- and Tlsuper is the corresponding super- 
Alfvenic Bernoulli pressure. Integration of these equations 
leads to 



Ps = n,„i-i-no, 

Ps - -^super + Hi , 



(31) 
(32) 



where Ho and Hi are integration constants. It follows from 
Eqs. (EB and OU that 



super - 111 ^ Ho - h Psub 

\2^0 



(33) 



All the above-mentioned relations and considerations are also 
valid for the case M^ ^ - 1 = {[f,g]a,/3f > 1- Only 
Eqs. ( l26l l. ( l27l i , ( ISTT i. and iT3[ are not valid then, along with the 
inequality Eq. ( l28l l. However, what does remain valid is that 



the magnetic field of the underlying MHS equilibrium can be 
amplified or weakened: 



B = B, 



Bs 



yl^s 



(34) 



super 



1 



with domains where |B| can be larger, and domains where \B\ 
can be smaller than Bs ■ 

We want to focus on purely sub-Alfvenic flows. In this case 
the transformation equations from the static to the stationary 
fields can be written as 





Bs 


^VfxVg 


(25) 


Ps 


= Ps(f,g) 




V 


= 


(26) 






(27) 




We are inter 



B ^VaxVfS^ 



Bs 



V^^' 



p^Ps- 2Plvr , 

MaBs 



^//op(1-M2) 



(35) 

(36) 
(37) 



We are interested in the fact that an astrosphere is termi- 
nated by a boundary between two diff'erent magnetic fields. For 
the magnetic field, this boundary is a tangential discontinuity or 
an encounter of two magnetic fields with a large gradient across 
that boundary. This boundary, the astropause, can therefore be 
regarded as a current layer, so we need additional information. 
We can deduce that the electric current density is also trans- 
formed by 

yuoj = V X (Va X V/3) 

= V[a,yS]/,, X (V/ X Vg) + [a,yS]/,,V x (V/ x Vg) 
= V[a,yS]/,y xBs+^iQ {a,l3]fjs 



^<i sign [[Q',/?]f.g] M^ 
sign [[Q',/?]f.g] 



[ V/ (VMa ■ Vg) - Vg (VMa ■ V/)] 



V 



V X (V/ X Vg) , 



(38) 



1-Mj 



which implies that even in the case of a static equilibrium with 
vanishing current density, i.e. a potential field, one gets a sta- 
tionary equilibrium with a non- vanishing current. 



__ sign[[Q',/?]fg]MA 
A'oV(l-^l)' 



J = 



[Vf(VMA-Vg)-Vg(VMA-Vf)]. 



(39) 



In addition to the previous works by the mentioned authors 
(Gebhardt & Kiessling,1992 ; Petrie & Neukirch [T999] l. we 
have found that the flow and the current are strongly correlated 
by means of the Alfven Mach number 

3. Two-dimensional equilibria 

To get exact and analytical equilibria, we restrict our view to 
two-dimensional equilibria since only symmetric equilibria are 
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known in infinite domains, see e.g. Tsinganos (ll982l) . There 
are no analytical and exact 3D MHD equilibria known that 
are bounded and that extend throughout the whole 3D space. 
Therefore, we assume from now on that f - A, where A is a 
function of x and y, and that g - z, so that we get the GSE (see 
e.g. Grad & Rubin fT958l l: 



dPs 



(40) 



and a relation between the Alfven Mach number and the deriva- 
tive of a: 



Mi = 1- 



1 



a'iAy 



« K(A))^ = (^) = ' 



dAj I -Ml 



(41) 



The current density for the stationary equilibrium can then be 
expressed by 



Aff = -^iQJ, = V ■ (Va) = a" (A) |VAr + Q''(A)AA 



(42) 



The calculation of the current can also be derived from 
Eqs. (l38T l and ( [39l ) where we have transformed the current den- 
sity for general 3D equilibria following Ampere's law, setting 
f - A and g - z- This is reasonable, if a domain in the vicinity 
of the equatorial plane of the star is to be represented by the 
calculated equilibria. 



4. The pattern of the flow field and the magnetic 
field: potential fields as magnetohydrostatic 
equilibra 

We show in this section that the distribution of (vutual or real) 
magnetic neutral points, which are also stagnation points if we 
assume a finite Alfven Mach number, determines the global 
topological and geometrical structure of the astrosphere. In the 
linear case, there is a relation between the magnetic multipole 
moments and the neutral points that will be given later on. But 
also for a more complicated class of nonlinear solutions of the 
MHS equations, it is possible to find similar relations. Those 
will be discussed in a future paper 

There are several reasons for using potential fields as ori- 
gins for our mappings. They are much simpler to handle than 
nonlinear MHS fields in the framework of our solution tech- 
nique, where the nonlinearity of the MHD equations is handled 
by a nonlinear mapping technique. Hence, one only has to solve 
a linear partial differential equation, whereas the nonlinearity is 
hidden in a nonlinear algebraic equation. In addition, potential 
fields have a highly stable character as they have no free mag- 
netic energy, although it is not known if stability is conserved 
after the mapping onto a stationary equilibrium (see Petrie & 
Neukirch ll999l l. Also, the connection between the neutral point 
distribution of the magnetic field and the global structure of the 
magnetic and the velocity fields can clearly be seen. 

Another reason is that there should be at least one saddle 
point (the so-called X-point) in this counterflow configuration. 
Amongst (A^ -H 1) null points, at least one point must exist in the 
vicinity of the nose of the astropause (at the stagnation point) 



for which the eigenvalues of the Jacobian of the Unearized mag- 
netic field 13, 
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take the form 







With the Grad Shafranov equation 
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dA 
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(44) 



(45) 



where 7(A) is the current function, we can write the linearized 
magnetic field as 



dx dy dx^ 



d^A 
"dx^ 



d\A_ 
dxdy 



Bi, 

Bly 



(46) 



The properties of the eigenvalues of the Jacobian of the mag- 
netic field determine whether the neutral (or stagnation) point 
is a point with which one can define a separatrix curve (or 
surface, see Arnold 1992 Reitmann 1996). The separatrix is a 
border surface between two different regions of a flow, or of a 
magnetic field that separates a vector field in areas of different 
topological connections. To get information about the topolog- 
ical structure of the magnetic field, we calculate 
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d^A _ ^ _ cl\A \ 

dx dy dx^ 
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Therefore, we obtain the following eigenvalues: 



d^A 
dxdy 

d^A 
dxdy 



d^A 
5x2 
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d\A 

5x2 



52a1 
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5x2 




52a 52a 1 


5y2 dx 


2 



(49) 
(50) 



For two conjugate complex solutions in 2D, we obtain A - 
+Ai € (C, i.e. only two purely imaginary values exist. If only 
two real eigenvalues exist, there is a saddle point (also called 
X-point), which is necessary for the existence of a separatrix 



~ The global coordinates (X, Y) are replaced here by (x,y) = (X - 
Xs,Y - Yg), with (Xs , Y^ ) being the coordinates of the nuU point. 
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and which guarantees that there is a boundary surface between 
two distinct areas of the flow; i.e. an astropause exists. If A is a 
potential field, i.e. AA = J{A) = 0, then the null point is a sad- 
dle point, as can be seen from Eq. ( |49] l. If the eigenvalues at the 
null point are purely imaginary, then, depending on the abso- 
lute value, a centre (a so-called 0-point) exists with topological 
circles as fieldlines, or a so-called /ocMi with spiral fieldlines. 
The last case is not found for solenoidal vector fields, since the 
trace of the Jacobian matrix vanishes for them. 

4.1. Potential fields 

We now have to calculate the flow pattern, i.e. the magnetic 
field pattern. This can be done by calculating the most sim- 
ple and stable magnetic fields, namely potential fields. We con- 
struct solutions by using a 2D multipole representation in the 
form of a Laurent series, which enables us to find static equi- 
libria. For this general kind of a conformal mapping, we ex- 
clude the region around the singularity {x,y) - (0,0) with 
Q - yjx^ + y^ < R,s within the termination shock, represent- 
ing the inner part of the astrosphere. We define u :- x + iy. 
J[ and S are the complex magnetic flux function and the com- 
plex magnetic field. J[{u) and S(m) are holomorphic functions 
(at least nearly everywhere), with a real part 51(^1) - (p,,, and 
imaginaiy pan 3(^) = A, and 5l(S) = B^ and 3(S) = -By ac- 
cordingly. A is the magnetic flux function and 0„, is the scalar 
magnetic potential. Therefore, we can write 



^_dm _ d,p,„ dA _ 
au ox ox 

We use the following Ansatz for the magnetic field 



B^Bsoo + Yj '^/^""^ ■ 



/i=i 



to satisfy the asymptotical boundary conditions 



Um S = Bsoc , 

|m|— »oo 

so 



-^ - BsaoU -I- Co In M -I- y CyU ^ 



(51) 



(52) 



(53) 



(54) 



is valid. 

There is a similarity between the logarithmic term of the 
hydrodynamical problem of a circular flow and the radial or az- 
imuthal part of a magnetic potential field, so that we can write 



r >F. 

-^ - Bsa,u + - — In M -H > CyU 

l.TTl ^—i 



(55) 



The first term in the expansion is the homogenous part due to 
the asymptotical boundary condition, which survives the non- 
canonical transformation. The second term is the circulation or 
monopole part due to a line current. If F = Fq (coSjSq + / sinySo) 
is real, then one has a typical counterflow configuration for the 
hydrodynamcal circulation of a flow around a circular cylinder. 
Here Fq - j-iqIo, where /q is the line current. The third term 
is a line-dipole part. In the case of pure hydrodynamics, this 
part represents the radius R^ = \Ci\ of a flow around a circu- 
lar cylinder, if sinySi - 0. The fourth term is the quadrupole 
part. Furthermore, |Ci| = 2/2fl is the dipole moment, I2 the cur- 
rent, a the half distance of the antiparallel line currents, and the 
product /2fl = const, while /2 — » 00 and a ^ 0. 

Moving stars, together with their winds, can be regarded as 
obstacles in the stream of the ISM. Such counterflow configura- 
tions lead to the formation of separating surfaces. On these sur- 
faces, stagnation points must exist where the flow velocity van- 
ishes. Stream lines passing through these stagnation points are 
called separatrices because they separate stream lines of differ- 
ent topological connection. They represent the borderlines be- 
tween two different flows. To calculate the stagnation point^ 
one has to solve the equation S = 0. In the topological the- 
ory of fluids, this field structure is called an X-point structural 
However, application to the Parker field would imply the exis- 
tence of an additional singularity beyond the X-point, located at 
the origin of the magnetized wind plasma, i.e. at the location of 
the central star itself. To generate a field similar to the hydrody- 
namic Parker flow of the heliospheric flow field (Parker ll96H . 
it is necessary (i) that xsp = xn{< 0) is the magnetic neutral 
point (where the magnetic field vanishes and in our treatment 
also the velocity field), (i) that for g = -\J(x^ + y^) -^ a ra- 
dial field structure exists, and (iii) that the flow field converges 
asymptotically to a homogenous field for £> — > 00. The Parker 
flow field is only useful and valid far away from the origin of 
the magnetized wind plasma of the central star One reason is, 
of course, that the field strength towards the origin is diverging. 
It turns out that selection of a circulation with a nonvanishing 
imaginary part is necessary. This enables the stellar wind to es- 
cape from the region of the reverse shock, and it creates a radial 
field structure towards the origin. 

Note: In our case, the Alfven or the usual Mach number 
does not determine the geometrical shape of the astropause, in 
contrast to the discussion in Parker (11 961) . Hence, our method 
cannot be compared directly with Parker's calculations. The 
mapping technique will allow sub- and superalfvenic flows to 
exist, although the streamline geometry does not change. Thus, 
it is also possible to have open field lines on both sides (upwind 
and downwind), like in the purely magnetic model of Parker 



v=l 



or especially 

ro(cos^o-H/sin^Q) , (cos^i+isin^ — 

^ — DsooM H — : in M — |Ui| -^ For a regular configuration, in the sense of finite Mach number 



-\C: 



2m 
(cosySi + isin/32) 



and density distribution, these nulls or magnetic neutral points are also 
stagnation points. 
+ terms of higher order . (56) 4 §^^.1^ p^^^jj, g^g gg^j^jjg p^j^jj, ^f ^j^g magnetic flux function. 
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4.2. Neutral points and multipole moments 

The complex magnetic field S can be calculated from ^ = 
(pi„ + iA, with (pm as the magnetic potential 



du 
where the Ansatz of the series ( l56l l yields 
r 1 Ci 2C2 



S = Bsoo + 



2ni u iP- 



(57) 



(58) 



the termination shock in the direction of the astrotail. We per- 
form a systematic calculation of the magnetic neutral points, 
making the general Ansatz of a Laurent series. If the Alfven 
Mach number is finite, the magnetic neutral points are identi- 
cal with the stagnation points. We analyse how the position of 
the stagnation points of the interstellar counterflow influences 
the field structure, especially the geometry of the astropause. 
We use the logarithmic part and the homogeneous asymptotic 
boundary condition, together with the assumption of a finite 
number of neutral points: 



A direct analytical method of calculating the null points 
should be applied here in the case of a multipole representation 
with two non vanishing multipoles. We restrict our analysis to 
the first two multipoles to get the magnetic null; 



-^ - BsoaU -H Co In M -H y CvU ^ . 



(64) 



r 1 Ci 

S = Bsoo + — ^=0. 

Zni u u^ 



This gives a quadratic equation, 
0, 



JT __Ci_ 

IjtBsoo Bsco 



(59) 



(60) 



Here v = 1 is the dipole, v - 2 the quadrupole, v = 3 the 
octopole, etc. The magnetic null or neutral points Uk are given 

by 



dJl 
du 



= 0. 



having the solutions 



Thus, we have to find the null points of the polynomial 
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v=l """ I<=1 

With the help of Vieta's theorem of roots, we get 

w+i 






+ 1 



■2 + J2 + 1R 
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r,r, 
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n «* 



1 < v<A^. 



(65) 



(66) 



(67) 



(68) 



(61) 



with 

Bsoo l6n^Bl^ 



J = 



Cu r,r, 

Bs^ in^BlJ 



(62) 



where the indices r and / indicate the real and imaginary parts 
of the coefficients. Therefore, 



Xs\ 



ys\ 



r, 



AnBs 



i(^ 



w +r-+'R 



The symbol CjJ'^j' denotes combinations of the Uk- These are the 
subset V -H 1 elements of the A^ + 1 elements (of the magnetic 
nulls). The distribution of the magnetic neutral points deter- 
mines the global geometry and the topology of an astrosphere. 

4.3. The equation of the astropause 

Regarding a symmetric and, with respect to the direction of 
interstellar medium, closed astrosphere, we only take neutral 
points on the x axis into account. The smallest, negative, x 
value gives 



_Jj_ 
AtiB, 



I 



Soo 



^2 V??2 + J2" + 2<R 



(63) 



A(Mi)=A,,pi =A(xi,0) 
and, therefore, 

A{x,y) ^Asep\u = 



^sep\ 



0, 



(69) 



(70) 



The singular point is situated where field lines meet; here is the 
beginning, i.e. the ending of several field lines. The stagnation 
point also marks one contour, e.g. a certain value of A. From 
that, we can calculate the equation of the astropause and the 
asymptotical equation of the astropause, which delivers the di- 
ameter of the astrotail at infinity. The second stagnation point, 
which is calculated for the case of the symmetric linedipole 
with ys2 - 0, xs2 > > jc^i, and xs2 > \xs\\, can be consid- 
ered as the radius of the inner astrosphere, i.e. as the position of 



for the separatrix with x <Q. However, including the monopole 
term, we obtain for the point where the separatrix (astropause) 
passes through the y axis: 



lim A{x,y) = A(x = -0,y = yo) = A.^pi 

A— >-0 



(71) 



where yo is the location of the astropause, which is lying on 
the y axis and, therefore, is elongated parallel to the inner astro- 
sphere; i.e. it is positioned at the same x coordinate of the star 
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Fig. 2. Field lines of the tail model as contour lines of the mag- Fig. 4. As in Fig.|2] but for U2 = 0.85 x 1 .5 = 1 .275. That some 
netic flux function and branches of the separatrix. ui is fixed at field lines do not appear closed is only a plotting artifact. 
-1.5, where the scale is in units of 100 AU. For U2 - 1.5 one 
can clearly see the similarity of this magnetic field lines with 
that of a flow around a cylinder 



> 




> 




Fig. 3. As in Fig.|2] but for U2 = 0.95 x 1.5 = 1.425. 

(xstar,ysrar) = (0, 0). Now, yo Can be determined from Eq. (fTTt . 
For jc > 0, we may calculate the curve of the astropause in the 
x-y plane as an implicit function of 



AH(x,y) = A(jc = -0,y = y£.)+7rCo 

= Asepl +7tCo - Asepl - 71 Co , 



(72) 



Fig. 5. As in Fig.|2l but for U2 = 0.5 x 1.5 = 0.75. That some 
field lines do not appear closed is only a plotting artifact. 

For X < and y > Q (with opposite sign for the limit of the 
arctan for y < 0, i.e. for the down part of the astropause), 

lim A{x,y) = A{0,yD) = Bsooyo - ^ Q, - — = = A,,,,{74) 
x^ -0 2 yo 

we get for the point where the separatrix intersects the y axis 

(75) 



r 



yo 



ttCo 



ttCo 



+ 



c, 



4-Bsco \\4-Ssoo/ Bsoo 
In the case of general asymmetric configurations, it is pos- 



taking the jump of the arcus tangens function into considera 

tion. Wifli three neutral points, the magnetic flux function can sible to calculate the intersection point with the y axis from 

be written as 

hmA(x,yD) = A(xi,yi) ^ A.^pi 
y Ixy 



3(^) = A = Bsooy+Co arctan - -Ci 



©-' 



■r 



1-C2- 



(x2 + y'^i 



.(73) 



Bs^yD - ^Cor + Co/ In lyol + ^ CvJd (76) 
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Fig. 6. As in Fig.|2] but for M2 = 0. In this special case the tail 
is similar to the Parker model (I196H . M2 = is not a null point, 
because this is the position of the singularity. M2 = is only due 
to the fact that in the sum of Eq. (l64l i only a monopole moment 
appears; see Eq. (l68T l. 




which is the part of the astropause that is opened in the tail di- 
rection. The outer separatrix is nothing else than the astropause. 
The asymptotical equation of the astropause is given by 

AHooix,y) := lim A(x,yH) = Bscoyn = ttCq => y// = ^ — -(80) 

In the series of Figs.|2]to|9]we show how the existence of two 
null points influences the shape of the astropause as an astro- 
spheric interface. The scale is in units of 100 AU. The first neu- 
tral point is chosen at a location that results in an astropause 
configuration applicable to the heliosphere. 

If we fix the neutral point, xi, in front of the astrosphere, 
we see that for X2 - -X] = R the fieldline geometry looks 
like that of a hydrodynamical counterflow of a cylinder with 
radius R (Fig.|2]l. By displacing the second null point towards 
the origin, i.e. for < X2 < -xi = R, the separatrix breaks up, 
and a tail-like channel is formed causing the drop-like shape 
of the astropause. This tail opens up for a null point approach- 
ing the origin (Figs. [3] - |5]l. For X2 = 0, the astrosphere has 
a typical Parker shape (Fig.|6). Further displacement to neg- 
ative values of X2 results in the formation of an anti-tailward 
bubble within the actual astropause (Fig.|7]i, leading to a notch 
in the astropause nose. This bubble grows until it touches the 
first neutral point resulting in a neutral point of second order 
(Fig. [SV. Here, the absolute value of the monopole moment is 
highesl3- A completely different scenario is shown in Fig.|9] 
Here, we have only one null point off' the x-axis. This implies a 
complex circulation F, leading to a spiral structure in the vicin- 
ity of the origin and emulating the Parker spiral. This represents 
an inner boundary condition giving a strong azimuthal compo- 
nent (a strong winding) of the magnetic field in the equatorial 
plane of the heliosphere. 

4.4. Discussion of tlie mirror symmetric case 

In the case of F, - and Cu = 0, we get an equilibrium that is 
mirror symmetric with respect to the x axis. The magnetic null 
points read 



-F, 



4nBs 



1^ 

as CO 



16n^ Bl 



Fig. 7. As in Fig.H but for M2 = -0.5 x 1 .5 = -0.75. Here, one One can easily see that, if 
gets a region with closed field lines in the downstream direc- _2 
tion. -£lL + i — 



<0 



Clr< 



F? 



16;r2 5 



(81) 



(82) 



With 

hm A(x,yo) = Bsocyo+Q) 

-v^+o 2 yo 



two magnetic nulls occur, which are not lying on the x axis, 
and the astrosphere is open with respect to the counterflow di- 

^sepi +■'!■ Co - ttCq, (77) rection. 

For 



we can calculate the other branch of the separatrix: 



Cu- = - 



A//(x,y) = BscoyD + ::Cor + Coiln\yD\+2_jCvyD 

^ v=l 

=" Bscoy + Co arctan ( - 1 - Ci -r. r - C2 

\xl x^ +y^ 

- ttCo, 



(78)- 



167^2 Bso 



(83) 



2xy 



Also that of the dipole moment. The fieldlines of that dipole es- 
cape to the right if the second null point is shifted from the right to the 

- left of the origin, so that the radial 'outflow' is stronger compared to 
/ 2 _i_ 2^^ tj I 

[^x +y ) tfje other images. The displacement of the second neutral point influ- 

(79) ences the inner boundary conditions for the field lines. 
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Fig. 8. As in Fig.|2l but for M2 
point of second order. 



-1.5; here we have a neutral 




Fig. 9. Field lines of the tail model as contour lines of the mag- 
netic flux function and separatrices for ui - -1.0 - i - (x,y) = 
(-1, -1); this is a potential field with a real valued circulation 
that shows the existence of a spiral-shaped field structure in the 
inner heliosphere, i.e. an azimuthal component; the thus im- 
itated Parker spiral should extend into the outer heliosphere. 
Due to a plotting artifact, not all calculated field lines are com- 
plete. 



only one stagnation point exists. The second null, m = 0, is not 
a real null point, but a pole. The first stagnation point is far 
away from the origin (i.e. from the star). The dipole part of the 
outflow, or the magnetic field, is like the resistance of a flow 



around an obstacle. Under the assumption Cir 
get 



u = X = 



AnB, 



+ R ^ 1 + 



r2 



Soo 



\6n^Bl R^ 

Soo 



BsooR\ we 



(84) 



For r, — > 0, we get two stagnation points (x - +R and x - -R), 
positioned on a circle with radius R. Figure |2] shows that the 
fieldlines of this magnetic field are identical to the image of 
streamlines of a flow around a body shaped like a circular cylin- 
der (keeping in mind the substitution 0,„ —>(/>, A —> i/r, with 
'W - (p + iif/ as the hydrodynamical potential and Bgco — » Vm). 
The circle is a separatrix that separates fieldlines of different 
topology. This astrosphere, however, would have the disadvan- 
tage of being a closed surface (line in 2D), and there would 
be no possibility that plasma could flow into the tail. This 
would imply a diff'usion process, as is described in Neutsch & 
Fahr (il983i) . 

5. Pure sub- or super-Alfvenic flows 

A tangential discontinuity has to form due to magnetic shear. 
Here two different, magnetized plasmas encounter At least, 
one has to expect a strong gradient perpendicular to the mag- 
netic field involving a non-singular current sheet. For typical 
astrophysical plasmas, the structure of such a current sheet 
can be derived by solving the coupled system of Vlasov and 
Maxwell equations selfconsistently, assuming the symmetry of 
the previous section. This kind of translation invariant plasma, 
although collisionless, can be considered to follow a quasi- 
Maxwellian distribution function, so that we can use the so- 
lution of Harris (11962) . Such a current sheet can be seen as 
the prototype of a current sheet separating two plasmas. To en- 
able a constant asymptotic homogenous field and to mimic two 
astropause current sheet positions, where the two symmetric 
branches of the potential field separatrices should be localized, 
we choose the following transformation equation 

A 



a{A) = 



sign [a'] A 



^ 



-1- fli In cosh 



Bs 



I -Ml 



-l-a2lncosh 
with derivative 



Bs, 



■yi 



d2 



(85) 



a'(A) = 



sign [a'] 

Jl - Ml 



ai 



Bscod 



■ tanh 



Bs. 



-yi 



SooUi 



02 



Bsood2 



tanh'^™ 



■yi 



d2 



fli 



B, 



S oo 



Bscod 



■ tanh 



A 



-yi 



fl2 



SooUi 



tanh'^^" 



di 



■yi 



, (86) 

DSooU2 (12 

where sign[ci;'] indicates that the asymptotical magnetic field 
can be parallel or anti-parallel to the asymptotical flow. Hence, 
with Eq. 



sign[a'] = sign[Boo]. 



(87) 



12 



D.H. Nickeler et al.: MHD flows in astrotails 



This transformation fulfills the symmetric asymptotic 
boundary conditions for the magnetic field, given by 



lim B - BooCv 



■■ lim B; 

J'— > — oo 



(88) 



where Boo is a constant. In the case of an MHS equilibrium, 
where 



A = Bsccy 



with 



Bsoc>Q, 



(89) 



the second term and the third term of Eq. (185) . which repre- 
sent the outer current sheets, are concentrated around yi and 
j2, which should be the locations of the astropause envelope 
borders in 2D (asymptotically). This implies that a' is the am- 
pUfication factor of the static magnetic field 



B - Va X e. = a'VA X e. = a' B, 



(90) 



which shows the behaviour of the asymptotic magnetic field 
if the magnetohydrostatic field is homogenous. The stationary 
field at infinity (x — > oo) can only depend on y because A con- 
verges to a value proportional to y and, therefore, a{A) ~ a{y). 
The symmetric boundary condition Eq. (l88l l leads then to 




Fig. 10. Regarding the field lines and especially the separatrix 
(thick line) one can see that they are not identical with the iso- 
contours of the electric current plotted in Fig. [in 



di d2 



(91) 



which means that the current sheets have different signs (an- 
tiparallel currents). The sign of the Jacobian (here a') must be 
unique in the whole domain, as there should be no roots for 
a purely sub-Alfvenic or super-Alfvenic flow. In addition, we 
make the assumption that < lii = d2 <^ \yi - yil and yi > y2, 
y\ = -yi > for the symmetric case (to ensure a highly sym- 
metric equilibrium). 



for Boo < « sign[a'] < 



2B, <B«(i+^1-M2^) 
2B, >B^(i-^1-M2^). 



(97) 
(98) 



5.1. Examples of noncanonical transformations for 
sub-Alfvenic flows 



With the above assumptions, Eqs. 
certain that 



The amplification factor a' of the static asymptotic magnetic 
field Bsoa at the axis of symmetry tells us that a magnetic jump 
occurs at the location of the separatrix, because the magnetic 
and Eq dgTT i we make fields from the inside and outside converge to different values. 
The strength of the jump is given by 



da 



dA 



= \a'\ > 1, 



(92) 



2Bi = Boo - a'iy = 0)B, 



5oo 



Ba 



SO that the flow is sub-Alfvenic. With condition ( |9TT l. Eq. ( |92] i 
can be written as 



V 



■Ml 



V^ 



m]. 

A,i J 



(99) 



\a'{y = 0)| 



sign [q-'] 



7 



-M] 



2Bi 
Bsoo 



Bco 

Bsoo 



IB I 
Bsoo 



> 1. 



Therefore, not all values of the inside magnetic field are al- 
lowed if the flow is purely sub-Alfvenic. The inside magnetic 
(9-3) field Bi is given by 



with 

Boo 

Bsoo 



sign [a'] 



Bi — Boo — 2Bi — Ba 



Vi 



m] 



11- Mj 



This leads to the following restrictions 
For Boo >0 « sign[Q''] > 



(94) 



7 



</ 



(100) 



Therefore, the polarity cannot change from the outside to the 
inside. With Eqs. dlOOb and ( |93] |. we get a lower limit for the 
inside magnetic field if the outside magnetic field Boo and the 
outside asymptotical Alfven Mach number M^oo are given: 



2B, <Boo(l-7l-^loo) 

(-7 



V 2Bi > B. 



Ml ] 

Aoo I 



(95) 
(96) 



|B,| > ^1 - Ml^ IBool = Bs« . (101) 

Boundary conditions for the ID case are given by the condi- 



tion Eq. ( llOOI l. This implies that the outer and inner Alfven 
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aurent density 




Fig. 11. Isocontours of the current density, where the high cur- 
rent isocontours are obviously concentrated around the he- 
hopause. The curvature and closure of the current isocontours 
can be seen clearly in contrast to the curved, but open separa- 
trix line in Fig. ( fTOl i: only in the asymptotical ID region, i.e. 
for ;ic ^ oo, the lines of maximum current and the separatrix 
(astropause) shapes seem to converge. 

Mach numbers are the boundary conditions for Eq. (IT4l . to- 
gether with the choice Bca = const. This determines B,, which 
cannot be prescribed as a boundary condition, since the prob- 
lem would be overdetermined. With the above relations we see 
that the possibility of setting boundary conditions is reduced, 
due to the reflection symmetry condition for the magnetic field. 
We choose the axis of symmetry in the x-y plane, the jc-axis as 
second boundary, and demand only regularity on the other two 
boundaries x = 1 and x - Xtaii-end- Regularity is guaranteed by 
the behaviour of the potential field and the transformation type 
in this domain. 

We can also use the above transformation to fulfill the 
boundary condition for the asymptotical ID region of a 2D 
field. In this case, the boundary conditions are mapped together 
with the mapping of the whole 2D potential field because we 
know in advance that this will again be a stationary equilib- 
rium state with field-aligned incompressible flow. Therefore, if 
for X — > oo, the equilibrium converges asymptotically to the 
ID equilibria given by Eqs. (fl4l l and (|86] |. this method can be 
used. We can take any of the given potential fields in the fore- 
going section, if we want to keep the potential character of the 
magnetic field in the tail. This leads then to a 2D sub-Alfvenic 
equilibrium state, writing 



]ima(A(x,y)) - a(Acc), with Acc-Bsooy- 



(102) 



In Figs. [To] and [TT] where we plotted bith the field lines and 
separatrix and the isocontours of the current density, an inter- 
esting feature of the transformation can be seen: while in mag- 
netohydrostatics, where B ■ VP = and j ■ VP = imply that 
the current is constant on field lines, so that the current isocon- 
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Fig. 12. Shown is the current density and its increase towards 
the inner astrosphere, i.e. the downwind region of the reverse 
shock. Also visible are current sheets with a width of 100 AU. 



tours coincide with the magnetic field lines, the situation is now 
completely difli'erent. Comparing Fig. [TT] with the fieldlines in 
Fig. [To] it can be seen that, asymptotically, the field lines and 
isocontours of the current density geometrically converge, but 
topologically they are different: while the field lines are open 
throughout the tail, the isolines of the current are closing in the 
vicinity of the separatrix (i.e. the astropause). 

To clearly show these differences in the field lines (Fig.[TO]i 
and the isocontours of the current density (Fig.fTTt. a finite 
width of the current sheet of lOOAU has been used. These 
broad current sheets are visible in Fig.[T2]where we have plot- 
ted the strength of the current density, whose absolute value 
increases towards the termination shock region. 

We now turn to the presentation of an example of a toy 
model that might represent the heliospheric tail region. First, 
we reduce the width of the current sheet to a more realis- 
tic value. That thickness can be estimated from the fact that 
it should be larger than several ion gyroradii (Neutsch & 
Fahr 1983; Fahr et al. 1982) . The ion gyroradius is at least of 
the order of 10^-10^ km, and we set the width of the current 
sheet to 10 AU. 

In addition, we can fall back on the measurements and es- 
timated values given in Frisch et al. (|E004) (Table 1, Model 
2) 15 For the outside magnetic field we take Boo = 5 //G, for 
the proton density n, ^ tie = 0.1 cm"^ and for the velocity 
Voo - 25 km s ' , so we obtain an interstellar Alfven Mach num- 
ber Maco ~ 0.72. Assuming an inner Alfven Mach number of 
about 0.52, the inner magnetic field in the vicinity of the x-axis 
becomes about 4juG. Taking the relation for the inner magnetic 
field, Eq. dlOOl ) and di = lOAU, we are able to calculate the 



^ The extreme value of the magnetic field is considered too high by 
P. Frisch, but preferred by Cox & Helenius (2003 1. 



14 



D.H. Nickeler et al.: MHD flows in astrotails 



X- component of the magnetic field 



Total strength of the magnetic field 





12.5 



12.5 



Fig. 13. The jc-component of the magnetic field. It rises towards 
the location of the termination shock caused by the magnetic 
monopole. Clearly visible are the jumps at the locations of the 
current sheets. 

y- component of the magnetic field 




Fig. 15. Tjtal strength of the magnetic field. The dominant con- 
tribution is from the x-component while, the contribution of the 
y-component is almost negligible. 

cmrent density 
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Fig. 14. The y-component of the magnetic field also shows a 
strong increase in its absolute value towards the termination 
shock region. 



12.5 



Fig. 16. The strength of the current density increases steeply 
towards the current sheets that have a width of 10 AU. The 
wavy shape of the peak current density is not physical but an 
artefact due to numerics. The current is normalized to units of 
2.65 X 10-1^ Am-2. 



transformation, using Eqs. (l86T l and (l9ll : 



where for Bi we used the definition given by Eq. ( |99] l: 



a'(A) 



Vi-< 



tanh 

Boo 



B, 4^^^^B. 



■y\ 



— tanh 



di 



^|^-MlooBo 



+^1 



Bi = i (Bo. - a'(y = 0)Bsc») = ^ 



\j Aoo 



(104) 



(103) 



The results for the difiisrent parameters are shown in 
Figs. fOUTTl in Figs.[T3] and [14] we plotted the x- and y- 
components of the magnetic field in units of 5 fiG. Towards the 
termination shock, the x-component grows especially around 
the jc-axis because it is approaching the mapped magnetic 
monopole. The jumps due to the current sheet are of the order 
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Fig. 17. This plot shows the Mach number and its strong gradi- 
ent across the hehopause. 

1 /iG. The y-component of the magnetic field is not symmet- 
ric with respect to the x-axis. Instead, a gradient arises due to 
the monopole. The contribution of the y-component to the total 
magnetic field strength is, however, very small, as can be seen 
in Fi g . [TSl compared with Fig.[T4li. Figure[T6]shows the strength 
of the current density. It is only different from zero at the loca- 
tions of the current sheets where it shows steep gradients. In 
the last figure (Fig.fTTIi we display the behaviour of the Alfven 
Mach number in the tail. It shows a strong gradient at the loca- 
tions of the current sheets. This gradient even increases towards 
the termination shock. 

6. Discussion and conclusions 

We present a method for calculating nonlinear MHD equilib- 
ria with an incompressible field-aligned flow. This method is 
applied to the scenario of a flow of interstellar plasma around 
the plasma bubble of a strong magnetized stellar wind. We use 
the classical method of conformal mapping of flows around an 
obstacle as the starting point of our calculations. 

We exclude violent structures (shocks) of flows in order to 
concentrate on the study of (i) the geometry of the contact sur- 
face, (ii) the surface currents that are coupled to the inner and 
outer magnetic fields and the Alfven Mach number as boundary 
conditions. The advantage of such a method is its high flexibil- 
ity in modelling the tail of stellar winds and the surrounding 
interstellar medium wind. What we need at least is informa- 
tion on the singular points (stagnation- and magnetic neutral 
points), their numbers, and their orders. Hence, it would be bet- 
ter if in situ measurement of the magnetic field structure could 
be made. Within the next decade, this is only possible for our 
own astrosphere, the heliosphere. 

In this paper, we restricted ourselves to thin nonsingular 
current sheets that have the special shape of a Harris-sheet or 
z-pinch configuration. The validity can only be justified within 
a multi-fluid theory or, better, within the framework of kinetic 
plasma theory and a detailed knowledge of the plasma envi- 



ronment in astrotails. Again, this aspect will be studied best 
observationally (in the next decade) for our heliosphere. 

As an improvement, additional current sheet structures 
should be taken into account, as e.g. the heliospheric current 
sheet is believed to extend beyond the heliospheric termination 
shock (see e.g. Pogorelov et al. . 120041 1. 

Our future aims are to use nonlinear static MHD equilib- 
ria as original equilibria, where the correlation between mag- 
netic neutral points will be much more complex. In addition, 
we have to find mappings with corresponding boundary condi- 
tions. Another important point is the symmetry we have taken 
into account: symmetry of such configurations can be broken 
easily by an angle between the magnetic field and the proba- 
ble flow direction in the vicinity of the heliosphere (see e.g. 
Frisch 1993). Such symmetry breaking will probably exclude 
an axially symmetric treatment of the problem. Nevertheless 
axially symmetric static equilibria can, in priciple, be used. In 
the case of pole-on counterflows of magnetized stellar winds, 
axial symmetry should be applied. 

For further investigation, transformations should be used 
that allow for transitions from sub- to super- Alfvenic flows per- 
pendicular to the magnetic field lines, as briefly described in 
Gebhardt & Kiessling (1 19921 1. Their application is much more 
complicated, as it is necessary to introduce four Euler poten- 
tials for the representation of the velocity and the magnetic 
field. 

Acknowledgements. This project was partly supported by the grant of 
the project No. 98030 of the PECS programme. D.H.N, thanks M. 
Karlicky for helpful discussions. 



References 

Arnold, V. I. 1992, Ordinary difi^erential equations (Berlin: Springer) 
Baranov, V. B., & Krasnobaev, K. V. 1971 Cosmic Research, 9, 620 
D'haeseleer, W. D., Hitchon, W. N. G., Callen, J. D., & Shohet, J. 

L. 1990, Flux coordinates and magnetic field structure (Berlin: 

Springer), 100 
Chandrasekhar, S. 1956, ApJ, 124, 232 
Cox, D. P., & Helenius, L. 2003, ApJ, 583, 205 
Fahr, H.-J., & Neutsch, W. 1983, A&A, 118, 57 
Fahr, H.-J., Grezedzielski, S., & Ratkiewicz, R. 1988, Annates 

Geophysicae, 6, 337 
Fahr, H.-J., Fichtner, H., & Scherer, K. 1993, A&A, 277, 249 
Fahr, H.-J., Neutsch, W., & Lay, G. 1982, Mitteilungen der 

Astronomischen Gesellschaft, 57, 294 
Frisch, P C. 1993, ApJ, 407, 198 
Frisch, P C. 2004, Adv. in Space Res., 34, 20 
Gebhardt, U., & Kiessling, M. 1992, Phys. of Fluids B, 4, 1689 
Goedbloed, J. R, & Lifschitz, A. 1997, Phys. of Plasmas, 4, 3544 
Grad, H., & Rubin, H. 1958, Proc. 2nd Int. Conf. on Peaceful Uses of 

Atomic Energy, Vol. 31, (Geneva: United Nations), 190 
Harris, E. G. 1962, II Nuovo Cimento, 23, 1 15 
Hameiri, E. 1998, Phys. of Plasmas, 5, 3270 
Imai, I. 1960, Rev. of modem Phys., 32, 992 
Kruskal, M. D., & Kulsrud R. M. 1958, Phys. Fluids, 1, 265 
Linde, T. et al. 1998, J. Geophys. Res., 103 (A2), 1889 
Lovelace, R. V. E., Mehanian, C, Mobarry, C. M., & Sulkanen, M. E. 

1986, ApJS, 62, 1 
Nemey, S., Suess, S. T, & Schmahl, E. J. 1991, A&A, 250, 556 
Neutsch, W., & Fahr, H.-J. 1983, MNRAS, 202, 735 



16 D.H. Nickeler et al.: MHD flows in astrotails 

Parker, E. N. 1961, ApJ, 134, 20 

Petrie, G. J. D., & Neukirch, Th. 1999, Geophys. Astrophys. Fluid 

Dynamics, 91, 269 
Pogorelov, N. V., Zank, G. P, & Ogino, T. 2004, ApJ, 614, 1007 
Reitmann, V. 1996, Regulare und chaotische Dynamik (Leipzig: 

Teubner) 
Scherer, K., Fahr, H.-J., & Ratkiewicz, R. 1994, A&A, 287, 219 
Schindler, K. 1972, in Earth's Magnetospheric Processes, ed. B. M. 

McCormac (Dordrecht: Reidel), 200 
Schindler, K. 1979, Space Sci. Rev., 23, 365 
Smith, E. J. 2001, J. Geophys. Res., 106 (A8), 15819 
Suess, S., & Nemey, S. T. 1992, in Solar wind seven, ed. E. Marsch & 

R. Schwenn (Pergamon Press), 281 
Suess, S. T, c& Nemey, S. T. 1990, J. Geophys. Res., 95, 6403 
Tsinganos, K. C. 1981, ApJ, 245, 764 
Tsinganos, K. C. 1982, ApJ, 259, 832 
Webb, G. M., Brio, M., & Zank, G. R 1994, J. Plasma Physics, 52, 

141 



D.H. Nickeler et al.: MHD flows in astrotails, Online Material p 1 



Online Material 



D.H. Nickeler et al.: MHD flows in astrotails. Online Material p 2 



Appendix A: Derivation of the transformation 
equations 

All solenoidal fields (i.e. vector fields with vanishing diver- 
gence) can be described locally by means of two scalar func- 
tions / and g 



[/' stfi V X (Va X V/J) X (Va X Vp) 
-i(VaxV^)2v[/,gt 



(A.9) 



B = V/xVg. (A.l) 

Applying the scalar product V/ or Vg, it follows that 

B-V/ = 0, and B-Vg = 0. (A.2) p^ = n s P +p|vp/2 



If we identify the equilibrium magnetic field with 

VaxV/3 = B, (A.IO) 

the sum of thermal and ram pressure with 



(A.ll) 



Since the potentials / and g are constant on field lines, 
besides the pressure in magnetohydrostatic equilibria Ps, we 
can understand f as a function of/ and g (i.e. intersections of 
/ = const and g = const are field lines). We now want to con- 
sider the magnetohydrostatic equations with Bs = V/ x Vg, 
where Bs is the magnetic field of a known static MHD equilib- 
rium. With the aforementioned equations, it follows that 



VPs = — (VxB)xB, 

and With VPs = ^V/+^Vg, 
df dg 

we can extract the equations of motion 
^ = Vg-Vx{VfxVg), 



dPs 
' dg 



V/-Vx(V/xVg) 



(A.3) 
(A.4) 

(A.5) 
(A.6) 



The scalar potentials are often called Euler potentials. Using 
unmatched Euler potentials a andyS, with / = f{a,[i) and g - 
g(a,p), it follows from adopting the functional determinant or 
Poisson brackets defined as 



[f,8]a.l3 = 



dfdg df, 



(A.7) 



da dp dp da ' 
that the following relation holds: 

Bs^VfxVg^ [/, g],ji VaxVp:^ [f, g^/, B . (A.8) 

Therefore, 

A^o VPs = (V X (V/ X Vg)) X (V/ X Vg) 

= (V X [/, g]a,fi Va X Vp) X [/, g]a^fi Va x Vp 

= (V[/, g]a,p X Va X V^) X [/, g],,fs (Va X V/3) 

+ ([/, ghfi^ X (Va X VP)) X [/, gU^ (Va X VP) 



= l^'^U^gVafi X (Vff X VP)\ X (Va X VP) 

+ lf 8]lj} V X (Va X Vp) X (Va x Vp) 
= [/, 8]l,fl V X (Va X V/3) x (Va x Vp) 

+ (^ V[/, g]lj, ■ iVa X VP)] (Vff X V/3) 



and the corresponding equilibrium current density of the sta- 
tionary equilibrium with 

Vx(VaxVp)=HQJ, (A. 12) 

we find 

Vn= l[/,g]2 C/xB)- -LB2v[/,g]^ (A. 13) 

We recognize the identical form of the last term on the right 
side of Eq. dA.lSI l with the right-hand side of the equation of 
motion (Eq. lA.14T ). which is 

Vn = 1(1 - MlXj xB)- -1^2 y^j _ ^2) ^ (^^^4) 



where we identify the Poisson brackets of Eq. (IA.13) with the 
Alfven Mach numer dependent expression in Eq. (IA.14l i. Then 
it follows that 



1-Ml^[f,g]i,>0. 



(A.15) 



For non-canonical transformations, the Poisson brackets 
has a non-constant value or a value different from unity. For 
canonical transformations, the Mach number is zero, or for a 
constant and non-unity value of the Poisson brackets, the Mach 
number turns out to be constant. It is therefore possible to map 
known solutions of the magnetohydrostatic equations with the 
help of non-canonical transformations into stationary solutions 
with a sub-Alfvenic flow. For Ma > 1, we find an analogy with 
the equations of incompressible stationary hydrodynamics (see 
Gebhardt & Kiessling |1992| l. Application of the same transfor- 
mation as in the case of the sub-Alfvenic equilibrium results in 
the same magnetic field and therefore in the same electric cur- 
rent. On the other hand, for Mach numbers with M^ > 2, we 
cannot find any sub-Alfvenic Mach number or solution. 

The pressure is going to be inverted with respect to the sub- 
Alfvenic pressure to become: 



n 



super-Alfvenic 



nw-n 



sub-Alfvenic ■ 



(A. 16) 



Here, Hh is a background pressure, which guarantees that the 
thermal pressure stays positive everywhere and fullfills the 
physical conditions (as e.g. vs > va), as well as the boundary 
conditions. 

As shown above, it is possible to map known solutions of 
the magnetohydrostatic equations via the non-canonical trans- 
formations into stationary solutions with sub-Alfvenic flow. If 
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there exists a non-canonical transformation / - f{a,/3) and with Aq as a constant of integration. The function A{x) can be 
g = g{a,/3) or a = a{f, g) and/3 = /?(/, g), then those stationary expressed in the framework of a 2D multipole expansion with 
fields are given by the help of the imaginary part of the line current, i.e. with the 

radial magnetic field (represented by the imaginary part of the 



B = Vff X V/? = —^ V/ X Vg = , Bs 

Ma (a if, g) ,/?(/, g)) ^ ^ MA(f,g) ^ 



^]l^i)Pia{f,g),P(f,g)) ^]tiop{f,g) 

Ma 



4 



Bs 



(l-M2)yLiop 



Ml , Ml 

n(/,g) - -r^lBp = Ps{f,g) r^ — \Bs 

2 ^ /^X 



P^p(f,g)- (A.17) 



Appendix B: Trajectories 

If / is a flux function, then / is constant on field lines. For 
AA - -yuo dP/dA - J(A) and with the help of the implicit 
function given by ;t' = x(^' j)' where xo = x(^' y) - const, it 
follows that A(x) = const. The equation x = x(^' y) - const 
therefore describes the bundle of the field lines of A(x)- 
The equation 

d^A dA 

^ = ttC^^)' + :r ^^ = --"0 ^^/^^ -J (AM) (B. 1) 

dx^ dx 

can therefore be regarded as the diff'erential equation for A as a 
function of ;^', if we have 

^ ., . = (B.2) 

in the case of a vanishing current function (i.e. a potential field), 
i.e. for 7(A) - 0. This results in a non-linear partial differential 
equation. 



6. 1 Radial magnetic fields 

For AA - -fio dP/dA = J(A) - and for radial trajectories, 
it follows that, for y/x - const on straight lines, the function 
A(x) - const. X is thereby regarded as a function of x and y, 
with xi^^y) '■- yl^- Application of the Laplace operator on 
A{x{x, y)) results in 

AA = A"(x)C- + -)+A'{x)^^Q 

xf* x^ x^ 

^ A"(x)(x^+l) + 2A'(x)X^O. (B.3) 

The second row of Eq. ( IB.3b can be expressed as 

^(A'(^)(l+/)) = 0. (B.4) 

Integration leads to 

A(x) - const ■ arctanOif) + Aq , (B.5) 



complex logarithm). 



